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ABSTRACT 

We report the discovery of a well-defined correlation between B-band face-on central optical depth 
due to dust, rjg, and the stellar mass surface density, /x*, of nearby (z < 0.13) spiral galaxies: log(r^) = 

1.12(±0.11) • log (^ MQkpc-^ ~ 8.6(±0.8). This relation was derived from a sample of spiral galaxies 

taken from the Galaxy and Mass Assembly (GAMA) survey, which were detected in the FIR/submm 
in the Herschel- ATLAS science demonstration phase field. Using a quantitative analysis of the NUV 
attenuation-inclination relation for complete samples of GAMA spirals categorized according to stellar 
mass surface density we demonstrate that this correlation can be used to statistically correct for dust 
attenuation purely on the basis of optical photometry and Sersic-profile morphological fits. Considered 
together with previously established empirical relationships of stellar mass to metallicity and gas mass, 
the near linearity and high constant of proportionality of the Tg — /i* relation disfavors a stellar origin 
for the bulk of refractory grains in spiral galaxies, instead being consistent with the existence of a 
ubiquitous and very rapid mechanism for the growth of dust in the ISM. We use the t b — /i* relation 
in conjunction with the radiation transfer model for spiral galaxies of Popescu & Tuffs (2011) to 
derive intrinsic scaling relations between specific star formation rate, stellar mass, and stellar surface 
density, in which attenuation of the UV light used for the measurement of star-formation rate is 
corrected on an object-to-object basis. A marked reduction in scatter in these relations is achieved 
which we demonstrate is due to correction of both the inclination-dependent and face-on components 
of attenuation. Our results are consistent with a general picture of spiral galaxies in which most of 
the submm emission originates from grains residing in translucent structures, exposed to UV in the 
diffuse interstellar radiation field. 

Subject headings: dust, extinction — galaxies: fundamental parameters — galaxies: ISM — galaxies: 
spiral 
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1. INTRODUCTION 

Broadband photometric imaging surveys in the 
UV/optical (e.g., SDSS (Abazajian et al. 2009), KiDS 
(de Jong et al. 2012), EUCLID (Laureijs et al. 2011), 
GALEX MIS) are, and will continue to be, a main source 
of information from which the physical properties of 
galaxies must be deduced. It is, however, a well known 
issue that the UV/optical emission of galaxies is strongly 
attenuated by dust and that this attenuation should 
be taken into account (e.g. Driver et al. 2007; Masters 
et al. 2010). This is particularly the case for late-type 
galaxies, which are usually much more gas- and dust-rich 
than early-type galaxies (as recently re-confirmed using 
Herschel-d&ta by e.g., Rowlands et al. 2011; Dariush et 
al. 2011; Bourne et al. 2012). Furthermore, detailed 
imaging studies of dust emission in the Milky Way and 
nearby spiral galaxies (e.g. Molinari et al. 2010; Bendo 
et al. 2011; Braine et al. 2010; Fritz et al. 2011) show 
that most of the dust is associated with large-scale 
structures in the neutral and molecular gas layers, which 
in turn causes the attenuation to depend heavily on disk 
inclination (Tuffs et al. 2004; Pierini et al. 2004; Driver 
et al. 2007). 

While it is generally agreed that the UV/optical 
emission of late-type galaxies must be corrected for 
dust attenuation, it has proven to be a challenge to 
measure the opacities of the disks, and various ap- 
proaches exist. The most powerful method is to utilize 
infrared measurements of dust emission in combination 
with UV/optical data, since the attenuating dust is 
heated by the UV/optical-emission it absorbs, and the 
bulk of this energy is re-radiated longwards of 60 /im 
in the far-infrared (FIR) and submm spectral range. 
Approaches utilizing this UV/optical-FIR/submm 
information range from semi-empirical ones, such as 
the IRX absorption estimator (Meurer, Heckmann & 
Calzetti 1999; Seibert et al. 2005; Johnson et al. 2007), 
via SED fitting using energy balance considerations (da 
Cunha, Chariot & Elbaz 2008; Noll et al. 2009; Serra 
et al. 2011), to radiation transfer modelling approaches, 
which explicitly calculate the UV/optical illumination of 
dust and the resulting FIR/submm-emission (Silva et al. 
1998; Popescu et al. 2000; Bianchi, Davies & Alton 2000; 
Gordon et al. 2001; Misselt et al. 2001; Bianchi 2008; 
Baes et al. 2010, 2011; MacLachlan et al. 2011; Popescu 
et al. 2011, PT11 hereafter). Unfortunately, such a 
coverage of the full UV/optical-FIR/submm SED is 
seldom, or only incompletely available for the population 
of spiral galaxies (despite their significant dust opacities) 
due to the scarcity of wide and sufficiently deep FIR 
surveys. Consequently, these methods can often only 
be applied to more massive spiral galaxies and starbursts. 

In the absence of FIR data semi-empirical methods, 
based solely on UV/optical-data, such as the UV- 
spectral-slope (3 or the Balmer decrement are often 
applied (Cardelli, Clayton & Mathis 1989; Meurer, 
Heckmann & Calzetti 1999; Calzetti 2001; Kong et al. 
2004; Seibert et al. 2005; Wijesinghe et al. 2011). These, 
however, often depend on either multiple UV bands or 
spectroscopy, information which is also often unavailable. 



In this paper we present a correlation between the B- 
band face-on central dust-opacity Tg and the stellar mass 
surface density /i* of late-type, non-AGN galaxies, and 
demonstrate that this can be used in combination with 
the radiation transfer model of PT11 to statistically cor- 
rect samples of late-type galaxies for the attenuation due 
to dust using only broadband optical photometry. The 
basis for this analysis is the overlap between the UV- 
optical-NIR/spectroscopic Galaxy and Mass Assembly 
survey (GAMA, Driver et al. 2011), and the FIR/submm 
Herschel-ATLAS (H-ATLAS, Eales et al. 2010) survey. 
After briefly recapitulating some of the key concepts of 
the PT11 model (§2) and describing the data and sam- 
ples employed (§3), we extract direct values of t b from 
the UV/optical- FIR-data for a sub-sample of local Uni- 
verse late-type galaxies jointly measured by GAMA and 
H-ATLAS in the H-ATLAS science demonstration phase 
(SDP) field and use these to derive and calibrate the 
correlation between and described in §4. In 

§5 we show that the correlation correctly predicts the 
NUV attenuation-inclination relations for complete sub- 
samples of GAMA late-type galaxies characterized in 
terms of /z*, thus demonstrating its predictive power to 
correct samples of late-type galaxies without available 
FIR data for dust attenuation using only broadband op- 
tical photometry. We then discuss our findings in the 
context of the production/injection/survival of dust in 
spiral galaxies and the scaling relations between specific 
star formation rate, stellar mass, and stellar mass surface 
density in §6, and summarize our results in §7. We use 
AB magnitudes throughout this analysis, and adopt an 
H a = 70 km s" 1 Mpc -1 , M = 0.3, fi A = 0.7 cosmology 
(Spergel et al. 2003). 

2. THE RADIATION TRANSFER MODEL 

In this paper we quantitatively link the characteristics 
of the attenuation of starlight in spiral galaxies to the 
mass of dust in the galaxies as directly determined from 
the FIR/submm integrated photometry. This approach 
mandates assumptions about the spatial distribution 
of dust in the galaxies. Here we utilize the radiation 
transfer model of Popescu et al. (2011, PT11 hereafter), 
which is applicable to a wide range non-starburst, 
late- type galaxies. We refer the reader to PT11 as well 
as Popescu et al. (2000); Tuffs et al. (2004) for a detailed 
description of the model, its parameters, and the work 
done to test its performance. Here we only supply a 
brief summary and detail its application to our data. 
The model consists of a bulge with an old stellar 
population and two exponential disks describing the 
distribution of old and young stellar populations as well 
as of diffuse dust located in these disks. This diffuse 
dust component can be seen as representing dusty 
structures (such as Cirrus) with a substantial projected 
surface filling factor. Emission from a dustless bulge is 
parameterized through the inclusion of a bulge-to-disk 
ratio (B/D) to accommodate a range of geometries along 
the Hubble sequence. In addition the model includes 
a clumpy dust component with an embedded young 
stellar population representing star-forming regions. 
The fraction of UV emission escaping from the regions 
into the diffuse ISM is given by a factor 1 — F (fixed to 
F = 0.41 for this analysis, following PT11). This model, 
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specifically the relative scale lengths and scale heights of 
the stars and diffuse dust in the exponential disks, has 
been calibrated on and fixed to the reproducible trends 
found in the local edge-on spiral galaxies analyzed in 
the radiation transfer analysis of Xilouris et al. (1999). 
As such, the wavelength dependence of the scale lengths 
is also fixed. 

In the PT11 model the opacity of the disk at a given 
frequency and position can be expressed in terms of the 
central face-on optical depth of the combination of the 
two dust disks at a reference wavelength (PT11 use the 

B-band at 4430 A), t S b . The value of Tg can be expressed 
as: 

4 = . (1) 

' s,d,B 

where Md us t is the total mass of dust in the galaxy, 
r s ,d.B is the scale-length of the exponential disk in 
the B-band, and if is a constant combining the fixed 
large-scale geometry and the spectral emissivity of 
the Weingartner & Drainc (2001) model. For the 

purposes of the work presented here the value of Tg 
must be derived from observable properties, hence Eq. 1 
must be re-expressed in terms of observational quantities. 

With the geometry of the model fixed, we express 
the physical scale length of the disk at the reference 
wavelength, r Si d.B, using the corresponding angular size 
at the redshift at which the source is observed. We have 
chosen to determine this angular size in the r band, 
which is less affected by the effects of dust attenuation 
than shorter wavelengths, while being less affected by 
noise than longer passbands, in particular the NIR, 
which may also suffer from systematic uncertainties (cf. 
e.g., Taylor et al. (2011)). 

The mass of dust Md us t is determined from the 
FIR/submm observations available from H-ATLAS. This 
data extends longwards of 100 ^m, thus predominantly 
sampling the emission by cold dust in the galaxy and 
warranting the assumption that this range of the SED 
can be reasonably approximated by a modified Planck- 
ian Sv(v) ~ v" B(is, T) with (3 = 2 (i.e the dust emission 
is not heavily affected by a warm dust component and/or 
stochastic heating processes broadening the peak of the 
SED). This allows Eq. 1 to be re-expressed as: 



rf - 



= A 



(1+*) 



3-/3 



B((l + z)naso,T ) 9: 



(2) 



,d,r 



with A = 6.939 

10~ 13 arcsec 2 J Jy _1 s _1 Hz -1 m~ 2 stcr -1 , 6 s ^ dyr 
representing the r band angular size corresponding to 
the disk scale length in arcsec, SV^so) representing 
the observed mono-chromatic flux density at 250 /xm in 
Jy, and B(v,T) representing a Planckian with units of 

WHz -1 m~ 2 stcr -1 , with a restframe temperature To. 
T , ^(^250)7 an d 9 s ,d,r will be derived from measure- 
ments of spatially integrated FIR/submm SEDs and 
optical morphologies in § 4. The numerical value of A 
has been calibrated using the detailed radiation transfer 
analysis results of the Xilouris et al. (1999) galaxy sam- 



ple. A detailed derivation of Eqs. 1 and 2 and their link 
to the PT11 model together with a detailed description of 
the numerical calibration of A is provided in appendix A. 

The choice of using 250 ^m is motivated by the tradeoff 
between using a measurement as far in the FIR/submm 
as possible, thus dominated by thermal emission of cold 
dust, and the sensitivity of the available data as discussed 
in §3. The values of t b depend somewhat on the fitted 
restframe temperatures of the modified Planckian fits via 
the temperature dependence of the derived dust masses 
as shown in § 4.1; the typical uncertainty in the temper- 
ature of ~ 1 K corresponds to an uncertainty of the dust 
mass of ~ 15%. 

3. DATA & SAMPLES 

The GAM A survey, currently comprising 11 band 
UV-NIR photometry (FUV, NUV, u, g, r, i, z, Y, J, 
H, K), and 300 — 800nm optical spectroscopy covering 
« 144 deg 2 to r petr o < 19.4 mag with a spectroscopic 
completeness of >96 percent, provides a large, statisti- 
cally complete, sample of local Universe galaxies with 
homogeneous photometry, spectroscopy, and ancillary 
data. This paper makes particular use of the stellar 
masses provided by Taylor et al. (2011) and single 
Sersic-profile fits performed by Kelvin et al. (2012), 
providing morphological information. For a full de- 
scription of the GAM A survey see Driver et al. (2011), 
as well as Baldry et al. (2009) for details of the input 
catalogue, Robotham et al. (2010) for details of the 
high-completeness tiling scheme, and Hill et al. (2011) 
together with Andrae et al. (2013, in prep.) for details 
of the photometry By design GAMA's coverage in the 
UV-optical-NIR is complemented in the FIR by the 
H-ATLAS survey (Eales et al. 2010) using the PACS 
(Poglitsch et al. 2010) and SPIRE (Griffin et al. 2010) 
instruments on board the Herschel Space Observatory 
(Pilbratt et al. 2010). H-ATLAS achieves 5cr point 
source sensitivities of 132, 126, 32, 26, and 45 mJy 
in the 100 /im, 160 /im, 250 /im, 350 /im, and 500 /im 
channels respectively. The details of the SPIRE and 
PACS map-making process are described in Pascale et 
al. (2011) and Ibar et al. (2010), while the catalogues are 
described in Rigby et al. (2010). The FIR data used in 
this analysis are taken from the catalogue of H-ATLAS 
sources matched to optical SDSS/GAMA sources 
using a likelihood-ratio method (Smith et al. 2011a), 
with a required reliability of >80%. This catalogue 
is defined by the requirement of 5-cr detections at 250 ^m. 

In the context of this work it is critical to select nor- 
mal, i.e. non-starburst, non-AGN, late-type galaxies as 
reliably as possible, as the radiation transfer model em- 
ployed is designed to treat, and is calibrated on, such 
sources, and is therefore only accurate for these. We 
first select those GAMA sources which have NUV and 
r band detections, reliable spectroscopic redshifts, and 
which have not been identified as AGN following the pre- 
scription of Kewley et al. (2006) 27 , limiting ourselves to 
a local Universe sample with z < 0.13 to avoid evolu- 

27 Almost all these sources have available Scrsic parameters and 
stellar mass estimates. Those that do not are excluded from the 
sample 
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tionary effects 28 . From these 14998 galaxies we select a 
population of probable spiral galaxies in the NUV-r vs. 
r-band Sersic index n r plane (Driver et al. 2012; Kelvin 
et al. 2012) as those with 

(NUV-r) < -22.35 log(n r ) + 11.5, (3) 

resulting in a subset of 11236 galaxies which we will 
refer to as the OPTICAL sample. 

The separator defined by Eq. 3 corresponds to the blue 
line in Fig.l. It was determined from GAM A sources 
with debiased visual morphological classifications from 
GALAXY ZOO DR1 (Lintott et al. 2008; Bamford 
et al. 2009), and is designed to maximize the number 
of GALAXY ZOO spirals correctly classified while 
simultaneously minimizing the contamination of the 
spiral sample due to GALAXY ZOO ellipticals based on 
a figure of merit maximization approach . As such, it is 
a conservative classification and will exclude some real 
spirals. We find a contamination of the resulting spiral 
sample by reliably classified ellipticals of 5-10% . 

From OPTICAL we select the subset of sources with 
> 3cr detections in the 100/im and/or 160/im channel in 
addition to at least a 5 a detection in the 250/zm channel 
(i.e. from the H- ATLAS SDP field where both SPIRE 
and PACS catalogues are currently available). We then 
calculate values of t b for these 79 sources as described 

in §4 excluding two with t b > 30 (which we take to 
be indicative of starburst systems). Finally we exclude 
three sources which have been visually classified as 
early-type galaxies by Rowlands et al. (2011), resulting 
in a sample of 74 late-type galaxies with detections in at 
least two FIR bands, referred to as the OPTICAL+FIR 
sample. 

4. THE OPACITY - SURFACE DENSITY 
CORRELATION 

4.1. Deriving opacities 

For each of the galaxies in the OPTICAL+FIR 
sample the disk opacity was calculated from Eq. 2 using 
knowledge of T and S250 (derived from Herschel data) , 
and the r-band angular exponential disk scale 9 s ^, r 
(derived from 8 e , ss ,r the r-band single Sersic effective 
size in arcsec (i.e, the half-light radius) catalogued by 
Kelvin et al. 2012). 

To derive T from the Herschel data we fit isothermal 
modified Planckians (/3 = 2) to all available data points. 
The requirement of a detection at 160 fim or shortwards 
allows the spectral peak of the dust emission to be 
well constrained. We find a median value of 23.2 K 
for T . The value of T is almost independent of the 
wavelengths at which the Herschel data are obtained, as 
the median temperature of sources with only a 160 fim 
PACS detection is 22.58 K, while that of sources with a 
PACS detection only at 100 /jm is 23.58 K. The median 
temperature of sources with PACS detections at both 
100 /im and 160 ^m is 23.35 K. The difference in median 

28 This limitation in redshift also is conducive to our confidence 
in the derived values of the Sersic parameter n, as meaningful mor- 
phological fits become more difficult at larger redshifts. 



dust temperature of galaxies between the 100 /zm-only 
and the 160 /im-only sample corresponds to an uncer- 
tainty in the dust mass of « 15%. The requirement of 
a datapoint at 160 fim or shortward does not appear to 
induce a strong bias towards warmer sources, since the 
median temperature of our sample is consistent with 
the mean value of 22.7 ± 2.9 K for blue galaxies with 
3 x 10 9 < M* < 3 x 10 n M Q (roughly comparable to our 
sample) found by Bourne et al. (2012) using a stacking 
analysis of H- ATLAS data on optically selected galaxies. 

Overall, we believe that the isothermal model con- 
strained by SPIRE data at A > 250 ^m and a PACS 
data point at 100 and/or 160 /im represents a robust 
method of determining dust masses using minimal 
assumptions, due to the decrease in the uncertainty of 
both temperature and amplitude arising from a data 
point constraining the peak of the dust emission, and 
because the wavelength coverage (A > 100 /im) misses 
any significant emission arising from warm dust in 
SF regions or from stochastically-heated small grains 
in the diffuse ISM (e.g., Popescu et al. 2000; Gordon 
et al. 2001; Misselt et al. 2001; Galliano et al. 2003, 2005). 

Using Eq. 2 to determine t b requires knowledge of 
the angular size corresponding to the disk-scale length 
in the r band. The relation between the observable 
single Sersic effective size and the disk scale-length of a 
spiral galaxy, however, is non-trivially influenced by the 
relative fraction of emission from the bulge and the disk 
as well as by dust present in the galaxy, with the former 
causing the ratio between 9 e . ss . r and 9 Sy d,r to decline, 
while the latter tends to cause sizes to be overestimated, 
increasing the ratio. Pastrav et al. (in prep.) have inves- 
tigated the combined dependencies of the ratio between 
# e . ss . r and 9 s ,d, r on bulge-to-disk ratio, dust opacity, 
inclination and wavelength, and provide their results in 
tabulated form. In this work we have self-consistently 
determined t b for the OPTICAL+FIR sample using 
Eq. 2 and the results of Pastrav et al. as detailed in 
appendix B. In doing so, we have used a bulge-to-disk 
ratio of B/D = 0.33 found to be representative of the 
generally earlier type, more massive spirals (Graham 
& Worley 2008), such as those in the OPTICAL+FIR 
sample. We caution, that this use of an average value of 

B/D will introduce uncertainty into the derived values 
f 

of t b as shown in appendix B, and will revisit our 
results when reliable bulge+disk decompositions based 
on higher resolution imaging of these objects becomes 
available. 

4.2. Deriving stellar mass surface densities 

We compute the stellar mass surface density /z* us- 
ing the physical radius r e , ss , r corresponding to the single 
Sersic effective radius in arcsec provided by Kelvin et al. 
(2012) and the GAMA stellar masses M* from Taylor et 
al. (2011) as 

M * " 27rrf >flS>r ~ 2nD\(z)0l ss , r ' 1 > 

where D 2 A {z) is the angular diameter distance corre- 
sponding to the redshift z. We note that the stellar 
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0.1 1.0 10.0 

sersic index n 



Fig. 1. — NUV-r vs. n r for non-AGN GAMA galaxies (z < 0.13) with NUV and r band detections. The color-coded sub-sample 
represents the overlap with GALAXY ZOO DR1 sources with debiased morphological classifications (Lintott et al. 2008; Bamford et al. 
2009) with the color representing the probability that the source is a spiral. The blue line represents an automated morphological selection 
designed to reliably select late-type galaxies with minimal contamination, calibrated on the GALAXY ZOO data (see GR12 for details) . 
Overplotted as stars are the sources detected in the H-ATLAS SDP field with SPIRE and PACS, with those in OPTICAL+FIR located 
left of the divide. 



masses predicted by Taylor et al. incorporate a single 
fixed prediction of the reddening and attenuation due to 
dust derived from Calzetti et al. (2000). Thus, expected 
systematic variations in reddening and attenuation with 
inclination, disk opacity and bulge-to-disk ratio are not 
taken into account. However, as discussed by Taylor et 
al. (see also Fig. 12 of Driver et al. 2007) the resulting 
shifts in estimated stellar mass are much smaller than 
the individual effects on color and luminosity, and 
should not significantly affect the relation between disk 
opacity and stellar mass surface density derived in this 
paper. Taking this, and other effects into account, the 
typical uncertainty in the stellar mass estimated by 
Taylor et al. is ~ 0.1 dex. 

We also note that Taylor et al. (2011) make use of a 
Chabrier (2003) IMF and the Bruzual & Chariot (2003) 
stellar population library, and that hence, any system- 
atic variations due to the choice of IMF or the stellar 
population library are not taken into account. 

4.3. The relation between opacity and stellar mass 
surface density 

t b is plotted against /z* for the OPTICAL+FIR sam- 
ple in Fig. 2. The data points are shown as symbols ac- 
cording to their Sersic index, with the color correspond- 



ing to the NUV-r color. Using a linear regression analysis 
taking the uncertainties in both /z* and Tg into account 
we find a power-law correlation between the two with 
X 2 /N DO f = 0.97 (iV DO F = 72) as 

log(r|) = 1.12(±0.11)dog ( - f* _ 2 )-8.6(±0.8) (5) 

\M kpc / 

depicted by the dash-dotted line in Fig. 2. The errors 
represent the 1-cr uncertainties in the regression analysis. 
The correlation is tightest for sources with NUV-r colors 
of ~ 3 with a slight increase in scatter for bluer and 
redder colors. This increase is likely due in part to 
the assumed B/D ratio in the determination of t b as 
detailed in appendix B, but may also represent a larger 
range in opacities for bluer possibly more irregular 
system and redder systems which may appear red either 
due to dust, or because they are more passive systems. 
There is also evidence for such a population of passive 

spirals, i.e, spirals with low Tg high as presented by 
Rowlands et al. (2011). 

We find that 50% of the sample lie within 0.14 dex of 
the correlation ( A Sj0 .5log(r^) = 0.14), comparable to 

the median measurement error for log(r^). Thus it is 
possible that a large fraction of the visible scatter is due 
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Fig. 2. — B-band face-on central optical depth r B vs. stellar mass surface density /i* for OPTICAL+FIR galaxies. Symbols are coded 
according to n and NUV-r color (see figure). The dash-dotted line represents the best-fit (Eq. 5). The median uncertainties in r B and 
/i* are depicted at bottom left. The inset depicts the dust mass (derived from t/ using Eqs. 1 and 2) as a function of stellar mass. The 
dotted line represents a reference value with a slope of unity and an offset corresponding to Mj us t/Mt = 0.003. Median errors are depicted 
bottom right. 



to measurement uncertainties. 

4.3.1. Range of applicability and limitations of the relation 

In order to understand the range of applicability of 
the correlation shown in Fig. 2 and given by Eq. 5, as 
well as to identify possible biases caused by the use of 
a FIR-selected sample in deriving this result, we have 
overplotted the distribution of the OPTICAL+FIR 
sample in the /x» vs. M* plane on that of the OPTICAL 
sample in Fig. 3. The OPTICAL+FIR sample covers a 
range of 7.6 < log(/i*) < 9.0 in //* more or less uniformly 
and can be deemed applicable in this range. Fig. 3, 
however, also clearly shows that the OPTICAL+FIR 
sample is strongly biased towards more massive sources, 
as shown by the positions of the purple stars. This bias 
arises from the the fact that the OPTICAL+FIR is 
defined by the sensitivity of the Herschel instruments. 
In spite of this clear bias in stellar mass, however, the 
sample does contain sources which provide a tentative 
coverage of the entire stellar mass range corresponding 
to the range in \x* as seen in the OPTICAL sample. As 
discussed in § 5, this bias in stellar mass does not affect 
the applicability of Eq. 5 to large samples of galaxies as 
constituted by the OPTICAL sample. 



The range of 7.6 < log(/i«) < 9.0 for which Eq. 5 
may be deemed applicable only provides a complete 
sample of galaxy stellar mass above 1O 9 5 M0. Thus, if 
the correlation is to be applied to samples of galaxies 
which need to be complete in stellar mass, these should 
be accordingly limited until the correlation can be 
calibrated for lower values of /x* (corresponding also to 
lower values of ) . 

Fig. 3 also shows that [i* appears to be loosely corre- 
lated with M*, with higher mass galaxies having larger 
values of //* . This raises the question of whether the t b 
- //* correlation shown in Fig. 2 is actually a relation 
between t b and M*. Using a linear partial correlation 

analysis of X = log(M»), Y = log(/i«), and Z — log(r^) 
on the OPTICAL+FIR sample we obtain the partial 
correlation coefficients rxy,z — 0.271, rxz.y — —0.047, 
and ryz,x = 0.611, however, implying that /i* is indeed 
the dominant factor in determining t b . 

Furthermore, as Eq. 5 presents a relation between 
two properties which are both inversely proportional 
to an area, we must ask ourselves whether the result 
is actually a spurious correlation due to noise in the 
size measurements. Due to the moderate redshift limit 
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Fig. 3. — fi* as a function of M, for OPTICAL (gray) galaxies with isodcnsity contours. Ovcrplotted are the OPTICAL+FIR sources 
(circles) and those sources which fulfill the morphological selection but only have SPIRE detections (stars). Dash-dotted lines indicate the 
range in fit, for which the correlation has been calibrated. The median errors on both properties are shown at bottom right. 



of z < 0.13, however, the uncertainties on the size 
determination are much smaller than the range in sizes 
found for a given value of M„, showing that the spread 
in values of /U* is mainly intrinsic. 

Finally, we emphasize that the quantitative accuracy 
of the relation given by Eq. 5 depends on the applicabil- 
ity of the large-scale geometry of the exponential dust 
disks as calibrated in PT11 to the range of late- type 
galaxies with 7.6 < log(/U*) < 9.0. 

4.4. The dust mass-to-stellar mass relation 

It is clear that Tq is akin to a surface density and 
requires measurements of both a galaxy's stellar mass 
and size to facilitate it's estimation. Nevertheless, a 
major underlying physical driver for the result presented 
by Eq. 5 is a roughly linear correlation between the 
mass of dust and stars in late-type galaxies of the 
OPTICAL+FIR sample. This is shown in the inset of 
Fig. 2, where dust masses (derived from the values of 
t b using Eq. 1) are plotted against stellar masses from 
Taylor et al. (2011). The dotted line depicts a slope of 
unity with a dust to stellar mass fraction of 3 • 10~ 3 as 
a reference value. 

4.4.1. Comparison with other Herschel results 



Several previous works have provided data on the dust- 
to-stellar mass ratio for different samples of galaxies, al- 
lowing quantitative comparisons with our results. Skibba 
et al. (2011) present stellar and dust masses for the galax- 
ies in the Herschel KINGFISH survey (Kennicutt et al. 
2011). For a sample of spiral galaxies of type Sa and 
later with M* > 10 9 M Q , comparable to our sample, one 
finds an average dust-to-stellar mass ratio of —3.02 ± 0.5 
(derived from Table 1 of Skibba et al. (2011) ), com- 
parable within errors to our result. Furthermore the 
dust-to-stellar mass ratio inferred by our data is com- 
parable within errors to that found for spiral galaxies in 
the Herschel Reference Survey (HRS, Boselli et al. 2010) 
by Cortese et al. (2012), as shown in Fig. 5 of Cortese 
et al. (2012) for individual morphological types and in 
Fig. 9 of Smith et al. (2011b) for all spiral galaxies in the 
HRS. The agreement is particularly good for the dust-to- 
stellar mass ratios derived for earlier-type spirals which, 
on average, are more massive and are likely to be more 
directly comparable to our sample (Cortese et al. 2012, 
their Fig. 5). In addition Fig. 5 of Cortese et al. (2012) 
also shows that the dust-to-stellar mass ratio is nearly 
constant for galaxies with morphological type Sa and 
later, especially for galaxies with Hi-deficiencies gener- 
ally indicative of residing in environments comparable to 
our sample (i.e., not being members of massive clusters). 

Finally we find that the dust-to-stellar mass ratio of 
~ 3 • 10 ~ 3 , is also in general agreement with that derived 
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by Dunne et al. (2011) (~ 2 x 1CT 3 ) for low redshift 
galaxies using all H- ATLAS SDP field sources, although 
we note that these authors have employed a different 
calibration of FIR/submm dust emissivity. 

The roughly linear slope of the relation between dust 
and stellar mass in the OPTIC AL+FIR sample also 
generally agrees well with the data for late-type HRS 
galaxies plotted in Fig. 8 of Smith et al. (2011b). The 
data show a slope which is slightly sub-linear over a 
large range in stellar mass extending down to below 
10 9 M Q . At galaxy stellar masses above 10 9 ' 5 M Q (more 
similar to our OPTIC AL+FIR sample), however, the 
data exhibit a slope which is considerably closer to 
unity. The results of Bourne et al. (2012), who find a 
correlation between dust and stellar mass based on a 
stacking analysis of optically selected sources, display 
similar properties, with the relation between dust and 
stellar mass steepening with increasing stellar mass 
(their Fig. 16). We note that Bourne et al. (2012) do 
not apply a morphological classification, but rather 
categorize their sample into blue, green, and red bins 
according g — r color which is likely to place some of 
the dusty edge-on spirals included in our sample in the 
green or even red bin in which their data displays a 
steeper, more linear slope. 

5. TESTING THE PREDICTIVE POWER OF THE 
STELLAR MASS SURFACE DENSITY - 
OPACITY RELATION 

The t b - [i % relation as formulated in Eq. 5 is funda- 
mentally an empirical result, linking the total dust mass 
per directly measured area to a direct measure of the 
stellar mass per unit area. Thus, the relation is largely 
independent of detailed assumptions about the distribu- 
tion of dust within the galaxy disks. Nevertheless, since 
we have formulated the previous analysis in the terms 
of the t b parameter of the radiation transfer model of 
PT11, we are in a position to directly and independently 
test the physical consistency and predictive power of 
Eq. 5 using any observable effect which is predicted to 
be a function of t b by the PT11 model. To this end 
we analyze here two quantities which are dependent on 
the amount and distribution of dust in galaxies, using 
a much bigger (and largely disjunct in terms of FIR 
detections) sample than used for the calibration of the 
t b - fi* relation, i.e., the OPTICAL sample. These 
quantities are the inclination dependence of attenuation 
of NUV emission from late-type galaxies (considered 
in § 5.2) and the scatter about the well known scaling 
relation between the specific star formation rate 
and the stellar mass M* (considered in § 5.3). We will 
show that even though the r B - /i„ relation is calibrated 
on a very limited portion of the overall population of 
galaxies it is applicable to the general population of 
spiral galaxies with 7.6 < log(/x*) < 9.0. 

5.1. Deriving attenuation corrections 

The radiation transfer model presented in Tuffs et al. 
(2004) and PT11 allows the inclination-dependent atten- 
uation of a spiral galaxy to be calculated for known values 



of t b , disk inclination id, escape fraction F, and bulge- 
to-disk-ratio B/D using Eqs. 17 & 18 from Tuffs et al. 
(2004) and the model predictions of attenuation coeffi- 
cients tabulated in PT11 29 . For the UV the value of 
B/D is of negligible importance as the UV emission is 
almost entirely produced in the disk, even for early-type 
spirals, and we have assumed F = 0.41 throughout as 
calibrated in PT11. Values of i& and t b leading to an 
object-by-object estimate of attenuation are found as in 
the following. 

5.1.1. Deriving inclinations 

The OPTICAL galaxy inclinations are calculated 
from the r-band single Sersic fit axis-ratios of Kelvin 
et al. (2012) as id = arccos((6/a) ss ), where (b/a) ss 
represents the single Sersic axis-ratio in the r band. 
These inclinations are then corrected for the effects of 
finite disk-thickness as detailed in § 3 of Driver et al. 
(2007), with an assumed intrinsic ratio of scale- height 
to semi- major axis of 0.12. 

5.1.2. Calculating attenuation- corrected magnitudes 

Using the catalogued stellar masses (Taylor et al. 

2011) and the measured values of 9 e , r (Kelvin et al. 

2012) together with Eq. 5 we estimate the values of t b 
for the entire OPTICAL sample. 

We then correct the NUV absolute magnitudes using 
the radiation transfer model (Tuffs et al. 2004, utilizing 
the aforementioned tables of attenuation coefficients in 
PT11) together with the disk inclination id and t b . 

5.2. The attenuation-inclination relation in the NUV 

Previous work (e.g, Tuffs et al. 2004; Driver et al. 
2007; Masters et al. 2010) has predicted and shown that 
the attenuation of UV/optical-emission in spiral galaxies 
is a strong function of inclination, with this effect being 
particularly pronounced at shorter wavelengths, thus 
severely influencing, for example, UV-bascd tracers of 
star-formation. This attenuation-inclination relation 
implies that the median observed absolute magnitude of 
members of a given late-type galaxy population should 
increase as a function of inclination. 

Here, given measurements of the inclinations, we 
use the attenuation-inclination relation to test the 
predictive power and physical consistency of Eq. 5 
by calculating the intrinsic absolute NUV magnitudes 
MnuVj corrected for attenuation as detailed in § 5.1. 
On an object-by-object basis the values of Mnuv will 
display scatter, at the very least due to an intrinsic 
spread in the galaxies' physical quantities. However, the 
median of an optically selected sample should no longer 
display an inclination-dependence after correction, if the 
transfer of UV radiation in galaxies is adequately de- 
scribed by the PT11 model, and the t b - ^» correlation 
given by Eq. 5 is representative of the late-type galaxy 

29 The requisite data specifying attenuation as a function of 
different wavelengths are available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsarc.u-strasbg.fr/viz-bin/qcatrj/A-l-A/527/A109. 
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population as a whole (in the according range of /U*). 
In particular, given the bias towards massive galaxies 
in the OPTICAL+FIR, the applicability of Eq. 5 to 
the galaxy population as a whole is by no means obvious. 

Fig. 4 shows the distributions of corrected and uncor- 
rected absolute NUV magnitude Mnuv as a function 
of inclination given as 1 — cos(i<j), for two sub-samples 
of OPTICAL defined by distinct ranges of /j*, thus 
corresponding to very different mean values of Tg. The 
samples are drawn from the range of observed stellar 
mass surface density /i» for which Eq. 5 is applicable 
(7.6 < log(//*) < 9.0; see §4.3 and Fig. 3), and cover 
the complete range of available galaxy stellar masses. 
The ranges in it* have been chosen to ensure that 
the samples are not affected by biases due to noise 
scattering sources into or out of the range in /i* for 
which Eq. 5 has been calibrated. We find a median 
value of Tg int = 2.26 for the galaxy sample with 

7.8 < log(/i*) < 8.3 and a value of t|, = 7.04 for the 
sample with 8.3 < log(/x*) < 8.8. The median value for 

both samples combined is = 3.10, while that for the 
entire OPTICAL sample with 7.6 < log(/j*) < 9.0, the 

range for which Eq. 5 has been calibrated, is = 2.42. 

As can be clearly seen the uncorrected samples (red 
points in Fig. 4) display a clear inclination-dependent 
dimming of their magnitudes, with the median magni- 
tude increasing (i.e. dimming) from the face-on case 
(1 — cos(id) = 0) to the edge-on case (1 — cos(id) = 1). 
In both ranges of it* the attenuation-corrected values 
of Mnuv-, derived as described in §5.1.2, are shown 
in blue. The corrected values of Mnuv no longer 
display a dependence on inclination, indicating that the 
correlation found using the OPTICAL+FIR sample is 
consistent with the independent observable presented by 
the attenuation-inclination relation, and with the radia- 
tion transfer model of PT11. This is also consistent with 
the supposition that the bias towards massive/bright 
sources in the OPTICAL+FIR sample, discussed in 
§4.3 and shown in Fig. 3 30 does not affect the correla- 
tion's applicability to the much larger OPTICAL sample. 

These conclusions are reinforced on a quantitative 
level by the agreement between the observed median 
distribution of the uncorrected samples and the pre- 
dicted inclination dependence of a fiducial galaxy with 
Tg corresponding to the median of the sample, and 
Mnuv corresponding to the median of the corrected 
bin- wise median Mnuv: depicted by the dash-dotted 
line in Fig. 4. In addition to the predicted dependence 
of attenuation on inclination, the difference in gradient 
of the attenuation as a function of inclination predicted 
by PT11 for r^> = 2.26 (corresponding to the range of 

7.8 < log(/j*) < 8.3) and t^> = 7.04 (corresponding to 
the range of 8.3 < log(tt») < 8.8) is also shown in the 
data. 

30 The bias of the OPTICAL+FIR sample towards bright 
sources is also visible in Fig. 4 where the uncorrected (green) and 
corrected (gray) values of Mpjuy for the galaxies in the OPTI- 
CAL+FIR sample in the appropriate range in /i„ arc ovcrplottcd 
and predominantly lie at the bright edge of the distribution. 



A large uncertainty in quantitative interpretations of 
the attenuation-inclination relation such as these arises 
from the difficulty of correctly classifying edge-on sources 
due to their intrinsic thickness and bulge component. 
This may cause these sources to be shifted towards 
lower values of inclination or to be absent from the 
sample. Indeed there is a hint that at high inclinations 
the sample may be slightly biased against low mass 
galaxies and that dust-rich spirals in general may appear 
very red at these inclinations leading to a possible mis- 
classification as ellipticals and a bias against edge-on 
systems. Nevertheless we believe that our results are 
not affected by strong, inclination-dependent, selection 
effects, as the samples for both ranges of ii* are essen- 
tially flat in 1 — cos(id). Furthermore, the distribution 
of K-band absolute magnitudes Mk (which are almost 
free of of dust attenuation) show no inclination depen- 
dence indicative of the presence of a strong selection bias. 

Both sub-samples defined in zi* display considerable 
scatter in Mnuv (after correction for attenuation) , with 
the average inter-quartile range being 1.6, respectively 
1.7 magnitudes. This scatter is much larger than can be 
accounted for by the scatter in t' b shown in Fig. 2. The 
range of scatter in Mnuv attributable to the scatter in 
Tjg, approximated by A S! o.5log(r^) as quoted in §4.3 
can only account for a range of 0.68 respectively 0.9 
magnitudes for the ranges 7.8 < log(tt») < 8.3 and 
8.3 < log(/ix*) < 8.8 respectively, even in the edge-on 
case. Additionally the inter-quartile ranges do not 
display inclination-dependence, as would be expected if 
the scatter were predominantly due to object-by-object 
variations in the dust opacity. Thus, the sample scatter 
appears to be dominated by the intrinsic scatter in 
Mnuv- The histograms of Mnuv m Fig. 4 clearly 
show that the corrected sample is more peaked and 
symmetrical with respect to the uncorrected sample, 
and that the large shoulder at fainter NUV magnitudes, 
a product of the inclination dependence, is largely 
removed after correction. This is the case both for 
the optically- and FIR-selected samples, while the 
remaining breadth of the distribution (especially for the 
FIR sample) reinforces the conclusion that the scatter 
in Mnuv is intrinsic. The remaining tail extending to 
faint NUV magnitudes can most likely be attributed 
to passive spirals, as presented e.g. in Rowlands et 
al. (2011), and to contamination caused by early- type 
galaxies (« 5-10%). 

Overall, we conclude that the inclination-dependent 
dimming of UV emission from galaxies in the complete 
optical sample can indeed be predicted using the relation 
between /j* and t|, calibrated on the subset of sources 
detected in the FIR. The consistency of the correlation 
with the PT11 model also lends confidence to the suppo- 
sition that the considerable shift in median magnitude 
due to the inclination independent component of the at- 
tenuation (« 0.9 and ~ 1.8 magnitudes, as predicted for 
galaxies seen face-on in the ranges of 7.8 < log(/i*) < 8.3 
and 8.3 < log(/i*) < 8.8 respectively) is also correct, 
as this is self-consistently derived together with the 
inclination-dependent component. This is investigated 




Fig. 4. — Uncorrected (red circles) and corrected (blue squares) values of Mnuv vs - 1 — cos (*d) f° r two sub-samples defined in fi* as 
stated in the figure. The samples include all values of M» present in the relevant ranges of /i* of the OPTICAL sample. Sources are binned 
in equal numbers (250) with the bin-wise median Mjvc/v an d 1 — cos^) depicted by solid black circles connected by solid lines, and error 
bars indicating the quartile boundaries. OPTICAL+FIR-sources are overplotted in green(uncorrected) and gray (corrected). The black 
dash-dotted line traces the inclination dependence predicted by the PT11 radiation transfer model for a fiducial galaxy with sample-defined 
median Tg (2.26 resp. 7.04, see figure), and median intrinsic Mjvt/y, defined by the corrected sample. Median random errors are shown 
shown top left, however, we expect increasing systematic errors in the determination of disk inclination at higher inclinations (see text). 
The histograms show the collapsed distributions in Mjmuv f° r the OPTICAL sample before and after corrections for dust attenuation (red 
and blue histograms respectively, with upper ordinate) and for the OPTICAL+FIR sample also before and after correction (green and blue 
hatched histograms respectively, with lower ordinate). 



further in § 5.3. 

5.3. The scatter in the specific star formation rate vs. 
stellar mass relation 

Although we have shown that the - /i* relation 
in combination with PT11 is effective at predicting the 
inclination-dependent component of attenuation it is still 
important to gain a quantitative measure of the efficacy 
of this technique in predicting the face-on component of 
the attenuation, which is not so directly probed by the 
analysis of the attenuation- inclination relation in § 5.2. 
Here we seek to achieve this by utilizing a fundamental 
scaling relation between physical quantities derived 



from UV/optical emission properties of galaxies where 
the intrinsic scatter between the physical quantities is 
sufficiently small as to be exceeded by the scatter in the 
observed quantities induced by dust attenuation. 

A particularly convenient scaling relation for this 
analysis is the well-known relation between specific 
star- formation rate, ip*, and stellar mass, M*, since, 
when derived from NUV magnitudes, the values of SFR 
used to construct ip* will be strongly dependent on the 
efficacy of the method used to correct for attenuation, 
whereas, as shown by Taylor et al. (2011) and discussed 
in § 4, the values of M* are much less affected by dust. 
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Here we convert from de-attenuated values of M^uv 
to V* using the conversion given in Kennicutt (1998) 
scaled from a Salpeter (1955) IMF to a Chabrier (2003) 
IMF as in Treyer et al. (2007). We note that the exact 
choice of conversion factor has no bearing on the analysis. 

In Fig. 5 the 0>* vs - relation is plotted for the 

OPTICAL sample limited to 7.6 < log(/x*) < 9.0 and 
M. t > 10 9 - 5 M (following the range of applicability of 
the t b — /i» relation given in § 4). To differentiate be- 
tween the effects of the corrections for the face-on and 
inclination-dependent components of attenuation we plot 
the relation three times: with no attenuation corrections 
(top left); with attenuation corrections as detailed in 
§ 5.1 , but with all inclinations artificially set to the me- 
dian inclination of the sample (bottom left); and with 
the corresponding full inclination-dependent corrections 
(top right). The expected trend of decreasing ip* as a 
function of M* is seen in all three panels. Comparison 
of the top left and top right panels shows, however, that 
the application of the inclination dependent attenuation 
corrections, in addition to inducing a overall systematic 
shift by a factor of 0.6 dex in ■0*, reduces the scatter in 
the relation, from 0.63 dex in the uncorrected relation to 
0.43 dex in the corrected relation 31 

This suggests a substantial predictive power, both of 
the Tg — jU* relation and the PT11 model, since we 
have applied an object-specific and large multiplicative 
correction to the NUV fluxes (by factors ranging from 
2.5 to 6.3 interquartile with a median correction of 3.8), 
yet have nevertheless succeeded in markedly reducing 
the logarithmic scatter in the ip* vs. M* relation 32 . 
Furthermore, comparison of the scatter in the partially 
corrected relation in the bottom left panel (0.57 dex) 
with the 0.63 dex scatter in the uncorrected relation 
in the top left panel shows that the total reduction 
in scatter is due not only to the correction of the 
inclination-dependent component of the correction, but 
also due to the correction of the face-on component of 
the correction. This is a strong indication that the zero 
point of the attenuation corrections (i.e., the face-on 
attenuation predicted by the PT11 model, which is 
the major contributor to the total attenuation) cannot 
be strongly in error. If the face-on component of 
the attenuation would have been independent of the 
stellar mass surface density, the large range of predicted 
face-on optical depths at a fixed stellar mass would have 
increased the scatter, rather than have diminished it. 

We note that the intrinsic scatter of 0.43 dex (in- 
terquartile) in the corrected relation of Fig. 5 (top right 
panel) will have substantial contributions from random 
errors. Major sources of this random uncertainty proba- 
bly arise from measurement uncertainties in the determi- 

31 All measurements of scatter were calculated as the difference 
between the quartiles of the distribution in ?/)», averaged over 15 
equal sized bins in M*, and weighted by the number of galaxies in 
each bin. 

32 As shown from the analysis of the multivariate relation be- 
tween Tg, M t , and /i* in § 4 the spread in face-on optical depth at 
a given M* arises from the large spread in disk radii for galaxies 
of a given M t , in conjunction with the close-to-linear vs. M* 
correlation 



nation of disk surface areas as well as from the estimates 
of galaxy stellar mass (~ 0.1 dex). In addition, we re- 
call that the galaxy sample will be contaminated at the 
5-10 % level by mis-classified spheroids. Furthermore we 
cannot rule out that there is some significant intrinsic 
scatter in the t b — [i % relation which would also induce 
a component of scatter in the corrected ij)* ~ relation. 
All this suggests that the intrinsic scatter in the V* vs. 
relation for spiral galaxies may be very low. 

5.4. Implications for the distribution and optical 
properties of grains in galaxy disks 

The success of the t b — /i* relation in combination 
with the PT11 model in predicting both the face-on 
and inclination dependent component of the attenuation 
in spiral galaxies has implications both for the spatial 
distribution of grains in galaxian disks as well as for the 
optical properties of these grains. 

Firstly, the quantitative consistency between the mea- 
sured dust surface density and the inclination-dependent 
attenuation of stellar light in disk galaxies, as predicted 
by PT11, is consistent with most of the dust in disks 
being distributed in structures sufficiently large to have a 
substantial projected surface filling factor. Furthermore, 
recalling that the t b — /i* relation is calibrated using 
measurements of the total submm flux, i.e tracing the 
total mass of dust in galaxies, the reduction in scatter 
about the 0* - relation induced by the application 
of PT11 points qualitatively towards most of the mass 
of dust in spiral galaxies (as traced in the submm) 
being distributed in diffuse, translucent structures, with 
a large fraction of the grains being exposed to UV in 
the diffuse interstellar radiation field as assumed by the 
PT11 model. 

In order to make this statement more quantitative we 
have in Fig. 6 plotted the mean interquartile range in the 
0* - M* relation as a function of attenuation corrections 
based on an effective value of the dust opacity parame- 
terized by a multiplicative scalar value \ as X ' t b- Hi 
contrary to the model of PT11 in which > 85% of the 
total dust mass is diffusely distributed, a large fraction 
of the dust mass measured in the submm were contained 
in compact, highly self-shielded regions, not exposed to 
the diffuse interstellar UV radiation field, the minimum 
in scatter about the 0* - M* relation should be attained 
for a relatively small value of x (i-e. X << !)• Instead we 
find that the minimum scatter is attained for x = 1-06, 
but that a range of x ~ 0.95 — 1.1 is not significantly 
distinguishable. This implies that, consistent with the 
PT11 model a fraction of > 85% of the total dust mass 
is distributed in diffuse, translucent structures. A more 
detailed analysis of the dependence of scatter in scal- 
ing relations, though beyond the scope of this paper, 
could in principle be used to fine tune model assumptions 
about the fraction of dust in clumps which are heavily 
self-shielded to UV light in disks of spiral galaxies, and 
thereby further improve estimates of the absolute level 
of -0, in the relation as well as the intrinsic scatter of the 
physical quantities. 

An analysis of the type performed here also has the 
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Fig. 5. — Specific star formation rate ip„ as a function of stellar mass M» for a subsample of the OPTICAL sample with 7.6 < log(/^») < 9.0 
and M* > 10 9 ' 5 M©. The relation is shown before correction for attenuation by dust (top left panel), after the full inclination-dependent 
correction, described in § 5.1, using the PT11 model in conjunction with the Tg — ^t* relation (top right panel), and after a partial correction 
using the procedure of § 5.1 but artificially setting a uniform inclination with cos^j) = 0.5 for all galaxies (lower left panel). The sources 
are binned in 15 bins of equal size in M*, with the median depicted by a filled circle, and the bars showing the interquartile range. The 
scatter in the relation due to the scatter in the NUV is reduced from the uncorrected to the fully corrected case. The intrinsic values of 
j/>* arc shifted upwards w.r.t. the uncorrected values. The linear gray-scale shows the number density of sources at that position, with the 
same scale having been applied to all samples. The median values and interquartile ranges are shown together in the bottom right panel. 
The uncorrected values are depicted by stars and a dash-dotted line, the values corrected at a fixed inclination of cos(i) = 0.5 are shown as 
inverted triangles and a dashed line, and the fully corrected values are shown as circles and a solid line. The bin centers have been offset 
by 0.01 in log(M„) for improved legibility. 
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Fig. 6. — Weighted mean interquartile range of i/>* as a function 
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of M« derived for fractions \ ■ t b of r g sampled in steps of 0.01. 
The minimum value of 0.43 is attained for \ = 1.06, however, it is 
not significantly distinguishable from that of \ = 0.95, . . . , 1.1. 

potential to empirically constrain the ratio between 
grain emissivities in the submm and UV/optical range. 
This arises because, while the estimates of rt on which 

the Tg — jU* relation is based are directly proportional 
to the dust emission coefficient in the FIR/submm, 
the amplitude of the attenuation corrections depends 
upon the dust emission coefficient at (in this case) 
UV wavelengths. Specifically, the demonstrated ability 
to correct for the inclination-dependent and face-on 
components of attenuation using Eq. 5, which was 



derived and calibrated using the FIR, is consistent 
with the ratio of the UV/optical and submm grain 
absorption coefficients being as described in the model 
of Weingartner & Draine (2001). 

6. DISCUSSION 
6.1. Dust production in spiral galaxies 
As already noted in § 4, the almost linear relation 
between the opacity of a galaxy disk, t* b , and the surface 
density of stellar mass, /it*, appears to arise mainly from 
a roughly linear relation between the extrinsic quantities 
of total stellar mass, M#, and total dust mass, M4. 
Although only established directly for galaxies with 
available FIR/submm measurements, the applicability 
of the same M* — Md relation to a statistically much 
more complete population of optically selected spiral 
galaxies from the GAMA survey is consistent with our 
analysis of the attenuation-inclination relation and ip* - 
relation presented in § 5. 

The physical origin of this link between M* and 
Md is far from obvious. On the one hand, the stellar 
mass is dominated by old, low mass stars which formed 
early on in the ca. 10 Gyr lifetime of a typical spiral 
galaxy. On the other hand, the main known process of 
injection of dust grains into the ISM is the condensation 



Dust Opacity - Stellar Mass Surface Density 



13 



of metals in the atmospheres of AGB stars on timescales 
of - 2 ± 1 • 10 9 yr (Dwek & Scalo 1980; McKee 1989; 
Morgan & Edmunds 2003; Ferrarotti & Gail 2006; 
Zhukovska, Gail & Trieloff 2008; Gail et al. 2009; Jones 
& Nuth 2011) 33 much shorter than the ages of spiral 
galaxies. Moreover, detailed modelling of the life cycle 
of refractory grains in the solar neighborhood (e.g., 
Jones, Tielens & Hollenbach 1996; Jones & Nuth 2011) 
predict that grains in the Milky Way are destroyed 
by SNe shocks in the tenuous ISM on timescales of 
~ 10 8 yr, much shorter than the timescale for the 
injection of dust from AGB stars, requiring that almost 
all observed refractory dust in the diffuse ISM must 
have been (re-)formed in situ soon after its destruction. 
This picture of rapid destruction and formation in the 
ISM is, however, difficult to reconcile with key physical 
and chemical properties of pre-solar grains as found in 
meteorites, most notably the segregation into separate 
populations of silicate and carbonaceous grains with a 
high abundance of minerals similar to those known to 
be produced in stellar sources. As discussed in detail 
by Jones & Nuth (2011), one is consequently confronted 
with a conundrum: either the grain destruction rates 
in the ISM have been grossly overestimated, allowing 
most refractory grains in the ISM to have an origin in 
AGB stars, or, alternatively, an as yet unidentified but 
very efficient mechanism exists that can convert gaseous 
metals in a low temperature and low pressure ISM into 
solid particles with the observed physical and chemical 
characteristics of interstellar grains. 

The close to linear relation between Md and M* 
underlying the t|, — /i* relation naturally favors the 
existence of a mechanism for efficient growth of refrac- 
tory dust out of the ISM, since any such mechanism 
would tap into the full reservoir of metals in the ISM, 
which are related to the integrated star formation over 
the lifetime of a galaxy. By contrast, if interstellar dust 
were mainly composed of longer-lived grains injected 
by AGB stars on timescales of ~ 2 ± 1 • 10 9 yr, one 
would expect the total mass of dust to be approximately 
proportional to the star formation rate (SFR) multiplied 
by a residency time, at least for systems with ages 
larger than the average dust destruction timescale 34 . 
Since the residency time should decrease with increasing 
SFR (since the frequency of destructive SN shocks 

33 Refractory grains have been observed to form in the metal- 
rich ejecta of core-collapse supernovae, which, alongside type la 
supcrnovae, have also been postulated to be major sources of in- 
terstellar grains. However, with the possible exception of the rem- 
nant of SN1987A (Matsuura et al. 2011; Lakicevic et al. 2012), 
FIR/submm measurements of cold unshocked ejecta in the cen- 
tral regions of prototypical young supernova remnants (SNRs) have 
shown that the ratio of solid state to gas phase ejecta is modest in 
comparison to the grain-to-gas ratio in the ISM (e.g., Tuffs et al. 
1997; Green, Tuffs & Popescu 2004; Barlow et al. 2010; Gomez et 
al. 2011). Given that, to escape the SNR, the condensates must 
traverse the very shocks postulated to be the main sink for refrac- 
tory grains in the ISM, it seems unlikely that supernovae are major 
primary sources of interstellar grains. 

34 In the very early evolution of systems dust mass may increase 
in parallel with stellar mass simply due to continuos injection of 
dust into the ISM driven by star formation. Only after the the 
age of the system increases to more than the average dust lifetime 
can the mass of dust be expected to be proportional to the SFR 
multiplied by a residency time. 



should be proportional to the SFR), and the SFR per 
unit stellar mass is known to decrease as a function of 
stellar mass, a strongly sublinear dependence between 
Md and M* would be predicted, even if an increase in 
the dust yield as a function of metallicity is taken into 
account. As such, if the origin of dust grains in the 
ISM were predominantly stellar, one would expect an at 
most very weak dependence of dust mass on stellar mass. 

If the mechanism for growth of grains out of the 
ISM implied by the slope of the - /j* relation was 
sufficiently prompt and efficient it would ubiquitously 
lead to a high fraction (i.e, of the order unity), r\ of 
all refractory elements being condensed into grains in 
the ISM of all spiral galaxies (as also inferred by, e.g. 
Dwek 1998; Edmunds 2001, and Draine 2009). To test 
whether 77 really does assume a universally high value in 
the ISM of local Universe spiral galaxies, we can make 
use of the well-established empirical relations linking 
stellar mass with gas phase metallicity and gas mass 
for this galaxy population. Specifically, the product of 
these relations will yield a relation between total metal 
mass and M„, which, by multiplying the metal mass by 
a constant value for 77, will predict a relation between 
total dust mass Md and M*. 35 . This relation can then 
be compared with the observed relation. 

Fig. 7 shows the predicted relations between M* and 
A^dust as derived using the mass-metalicity relation 
for late-type galaxies (Tremonti et al. 2004; Kewley & 
Ellison 2008), converted to gas-phase metalicities and 
a Chabrier (2003) initial mass function (IMF) as in 
Peeples & Shankar (2011, PS11), and the gas-to-stellar 
mass ratio from PS11. The relations for 77 = 0.5 and 
77 = 1 are respectively shown by the solid black and 
dash-dotted gray lines, together with horizontally and 
vertically striped regions indicating the 1-a scatter 
around the relations. It is apparent that the observed 
trend in vs. Md ust , shown by the ovcrplotted data 
points from the OPTICAL+FIR galaxies, is indeed 
quite well predicted by the mass-metalicity and gas mass 
vs. stellar mass relations for constant 77, and that the 
required value of 77 indeed has to be high. If fact values 
of between 0.5 and 1 are required for 77, about a factor 
of two higher than the several tens of percent of ISM 
metals that are predicted to be present in the form of 
grains by several detailed physical models such as those 
by Dwek (1998, r\ ~ 0.4), Edmunds (2001, 77 > 0.4), and 
Galliano, Dwek & Chanial (2008, 77 ~ 0.3). However, 
as we discuss in Appendix A, it must be born in mind 
that the measured dust masses plotted in Fig. 7 were 
derived from the FIR/submm observations using a mass 
absorption coefficient which is actually quite uncertain. 
In particular, whereas the relative values of the UV and 
submm grain absorption cross sections of Weingartner 

35 The underlying assumption of such a prediction is that all 
galaxies will have experienced a similar star formation history. 
Variations in this history can give rise to significant scatter around 
the relation, in particular, the time at which a large burst of SF 
occurs may strongly influence the observed dust-to-stellar mass ra- 
tio. Nevertheless, as we have endeavored to select a pure sample of 
normal spiral galaxies, and the specific SFRs obtained for the OP- 
TICAL+FIR sample don't display bimodality, such an assumption 
does not appear unreasonable. 
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& Draine (2001) used by our radiation transfer analysis 
connecting the submm emission and UV attenuation 
characteristics of spiral galaxies have been empirically 
constrained with respect to hydrogen gas column 
through measurements of extinction and emission of 
diffuse dust in the Milky Way, the absolute value of the 
absorption cross section per unit grain mass /c m , needed 
to deduce the value of rj, is relatively uncertain. As 
noted by Draine et al. (2007) the value of n m for the 
model of Weingartner & Draine (2001) requires more 
heavy elements than appear to be available and the 
mass of dust may be overestimated by a factor of ~ 1.4. 
36 Such a shift in K m (while leaving the grain absorption 
cross sections relative to hydrogen unchanged) would 
reconcile the majority of the measurements plotted 
in Fig. 7 to an n of ~0.5, given the observed scatter. 
At the same time this would preserve the observed 
quantitative connection between the UV attenuation 
and the observed surface density of submm emission, as 
predicted by the PT11 radiation transfer model. 

The only way we could envisage maintaining this 
demonstrated ability to link attenuation of starlight to 
the observed surface density of submm emission, while 
avoiding rj approaching unity and avoiding having to 
raise K m for dust in the diffuse ISM from the values given 
by Weingartner & Draine (2001), would be to invoke 
a population of highly self-shielded compact dense 
clumps as the source for a large fraction of the observed 
submm emission from spiral galaxies, in conjunction 
with us having systematically overestimated the intrinsic 
angular sizes of the disks seen in r-band. 37 While the 
present accuracy of measurement of intrinsic disk sizes, 
as outlined in Appendix B, may not completely rule out 
such a scenario, recent high angular resolution submm 
imaging of the galactic plane of the Milky Way by the 
Herschel Space Observatory, sensitive to emission on all 
angular scales, clearly show that the vast majority of 
submm photons originate from translucent large scale 
structures (e.g. Molinari et al. 2010). 

We conclude that the near linearity and high constant 
of proportionality of the Tg — /x* relation, itself based 
on a near linear relation between M<j U st an d M*, is 
indeed in good agreement with a roughly constant and 
high (but still physical) fraction of all ISM metals being 

36 This is also reflected by the fact that determinations of re m 
based on metal abundance as an input constraint predict higher 
submm grain emissivities than those in the Weingartner & Draine 

(2001) dust model. For example the emissivity model used by 
Dunne et al. (2011) is partly based on the analysis of James et al. 

(2002) , who, by assuming that 45.6% of all metals are converted to 
dust, derived a dust absorption coefficient per unit mass at a wave- 
length of 850fim which is ~ 70% larger than that of Weingartner 
& Draine (2001). 

37 This would reduce the mass of dust needed to explain the ob- 
served submm fluxes, since the dust in self-shielded clumps, while 
no longer able to efficiently participate in the attenuation of light 
from stellar populations not spatially correlated with the clumps, 
might be expected to have a much higher value of re m in the submm, 
due to the formation of ice mantles in such environments. The cor- 
responding reduction in the mass of dust in the diffuse ISM would 
then need to be exactly compensated for by the reduction in the 
inferred intrinsic angular size of the disk, such as to restore the 
opacity of the disk to the levels needed to predict the attenuation 
of the starlight (as quantified through the attenuation-inclination 
relation and the scatter in the tp — M* relation). 
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Fig. 7. — Predicted values of dust mass M,j us t as a function of 
stellar mass M» for an assumed conversion of a fraction n of all 
ISM metals to dust. The mass-metallicity relation (Tremonti et 
al. 2004; Kewley & Ellison 2008) converted to gas-phase metal- 
licities, a Chabrier (2003) IMF as in PS11 and the stellar-to-gas 
mass ratio (PS11) used in deriving the expectations are overplot- 
ted as a dashed and dotted lines, respectively with the shaded 
areas indicating the range of 1-<t scatter around the relations. The 
predicted relation and l-cr scatter (derived as sum quadrature) be- 
tween Md us t and M* is shown for r\ = 0.5 ( solid black line and 
horizontally striped region) and for r\ = 1. (dash-dotted gray line 
and vertically striped region). The diffuse dust masses of the OP- 
TIC 'AL+FIR sample, derived from the values of r g using Eq. I arc 
overplotted as filled circles with error bars (errors on M<j uat take 
into account errors on Tg and f e ,ss,r) 

present in the form of grains, and is best understandable 
in terms of the existence of a ubiquitous and very 
rapid mechanism for the in situ growth of grains in the 
gaseous ISM. Based on a joint consideration of measured 
dependencies of dust mass, gas fraction and mctallicity 
on stellar mass, our simple analysis provides a direct 
and model-independent empirical confirmation of work 
which has used more sophisticated chemo-dynamical 
simulations of the dust cycle in local and high-z galaxies 
applied to dust abundances and gradients to infer a 
dominant in situ source of interstellar grains both in 
local, normal galaxies and high-z starbursts (e.g. Dwek 
1998; Calura, Pipino & Matteucci 2008; Draine 2009; 
Michalowski et al. 2010; Dunne et al. 2011; Dwek & 
Cherchneff 2011; Dwek et al. 2011; Inoue 2011; Matt- 
son, Anderson & Munkhammar 2012; Valiante et al. 
2011). Since we make no assumptions about the stellar 
populations other than that of the current injection rate 
of Stardust being proportional to the recent SF rate, 
our conclusion that Stardust is a minor constituent of 
dust in spiral galaxies holds even if the initial mass 
function for stars were to be top heavy, which has been 
suggested (e.g. by Dunne et al. 2011) as a possible way 
of alleviating the need for grain growth in the ISM. 

Moreover, the ability of the t b - /z* relation to predict 
the NUV attenuation-inclination relation suggests that 
the majority of grains are exposed to non-ionizing UV 
light in the diffuse interstellar radiation field, so are 
refractory in nature 38 and reside in the diffuse ISM. As a 
consequence, our results not only require a very efficient 

38 Volatiles in the form of ices will almost instantaneously return 
to the gas phase through photodesorption if exposed to UV in the 
diffuse interstellar radiation field 
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grain formation mechanism, but also require that this 
mechanism pertains to the formation of refractory 
grains, rather than merely to the condensation of ices in 
highly self-shielded regions. 

The nature of the mechanism for forming refractory 
grains in the ISM is completely open. In their compre- 
hensive analysis of the evolution of the interstellar dust 
population in the solar neighborhood in the Milky Way, 
based on a one zone chemical evolution model account- 
ing for the growth of individual species, Zhukovska, Gail 
& Trieloff (2008) conclude that the interstellar dust 
population is dominated by refractory grains grown by 
accretion of gas phase metals in dense molecular clouds, 
with Stardust (in their model from both from AGB 
stars and from supernovae) constituting only a minor 
fraction. This result is consistent with our conclusions, 
independently inferred from the t b — /i* relation for 
spiral galaxies, but only provided a mechanism exists 
to propagate the refractory grains from the clouds into 
the diffuse ISM on timescales shorter than the timescale 
for grain destruction in the diffuse ISM. Alternatively, 
Draine (2009) has proposed that refractory grains can 
grow in diffuse interstellar clouds, in the presence of UV 
radiation. This would seem to be more easy to reconcile 
with our result that the bulk of all grains must reside 
in translucent structures illuminated by UV, as it would 
bypass the need for a rapid propagation mechanism. 

In general, the r B — /i„ relation may be useful as a 
diagnostic tool to investigate the universality and nature 
of the in situ grain-formation mechanism. In particular, 
although our analysis favors grain condensation from the 
ISM as the main grain injection mechanism, our present 
statistics cannot rule out that a significant minority of 
the grains have a stellar origin. Jones & Nuth (2011) em- 
phasize that there is a considerable uncertainty in theo- 
retical predictions for grain lifetimes, so that, while there 
is a strong requirement for Carbonaceous grains to be 
rapidly recycled in the ISM, this requirement may be 
less strong for Silicate grains. Analysis of the t b — /i* 
relation for larger statistical samples will allow separate 
relations to be established for spiral galaxies as a function 
of recent SF history, spiral arm coverage (i.e., lateness) 
and specific SF rate which may throw more light on this 
question, particularly if accompanied with data on the 
strength of the 2200 A absorption and MIR Polycyclic 
Aromatic Hydrocarbon (PAH) emission features, both 
of which specifically probe Carbonaceous particles. 

6.2. The attenuation of starlight in spiral galaxies 

Having discussed the use of the t b — /!„ relation as 
a diagnostic of physical processes driving the efficient 
production of interstellar dust in the disks of spiral 
galaxies, we return to the main goal of this investigation, 
namely the use of the t b — [i % relation, in conjunction 
with the radiation transfer model of PT11, to correct 
for the attenuation of stellar light by dust. As described 
in § 5, this can be done on a object-to-object basis 
for large statistical samples of spiral galaxies, using 
readily available optical photometric properties for 
each galaxy. Although the predicted attenuations are 
quite substantial, especially in the UV, our quantitative 



analysis of the attenuation-inclination and of the scatter 
in the ip* - relation lends some confidence that the 
corrections are not largely in error. 

As we have already noted in the case of the V* - M* 
relation the scatter in fundamental scaling relations 
based on UV and optical quantities can be significantly 
reduced through application of attenuation corrections 
based on the t b — ^i* relation even when no dust 
emission data is available. The 0.43 dex interquartile 
scatter in the ip* - M* relation after correcting for at- 
tenuation already implies a very tight relation between 
the current and past star-formation in spiral galaxies in 
the local Universe that would need to be reproduced by 
any theory of the formation and growth of spiral galaxies. 

Although affecting the scatter of the V'* vs. M* 
relation, the attenuation corrections do not, as already 
noted, strongly affect the slope of this relation, at least in 
the range of M* for which the attenuation corrections are 
presently available through the t b — /i* relation. This is 
because, as shown in § 4, opacities are statistically much 
more tightly related to stellar mass surface density, 
rather than to stellar mass, coupled with the fact that 
spiral galaxies of a given stellar mass exhibit a wide 
range of disk sizes. However, this situation will no longer 
apply to scaling relations function of the stellar 

surface mass density, To illustrate this, we plot in 
Fig. 8 the relation between ip* and /z* , for the same sam- 
ple as used for the tp* vs. M« relation, both before and 
after correction. Remarkably, the slope in the relation 
between the uncorrected quantities is entirely removed 
after correction for dust, showing that ip* is statistically 
independent of /i„ . The scatter in the relation is reduced 
from 0.62 to 0.49 dex. The latter value is somewhat 
larger than that for the V* vs - relation, perhaps 

implying that the ip* vs. M* relation has a smaller in- 
trinsic scatter and is thus the more fundamental relation. 

Overall, the fact that disk opacities scale systemat- 
ically with stellar mass surface density, as opposed to 
being randomly distributed, may help to explain the 
preservation of systematic and in some cases surprisingly 
tight relations between optical or UV tracers of physical 
quantities, even when these observational tracers are 
heavily affected by dust attenuation, and may help 
to explain why many relations were historically dis- 
covered only with relatively crude corrections for dust 
attenuation. Apart from the scaling relations analyzed 
here, a further relation which would be particularly 
pertinent to reanalyze would be the Tully-Fisher (TF) 
relation between luminosity and dynamical mass (Tully 
& Fisher 1977), which is even tighter than the V* 
vs. M* relation. A similar analysis to that applied 
to the V* vs - relation here applied to the TF 

relation could provide a still sharper tool for statistical 
analysis of attenuation corrections, as well as potentially 
improving the accuracy of the TF relation at shorter 
wavelengths, both as a distance indicator, and as a di- 
agnostic of the formation and evolution of spiral galaxies. 

Finally, we re-emphasize that the corrected relations 
are only for subsets of the galaxy population restricted 



16 



Grootes et al. 



-9.0 — i — ■ — ' — ■ — i — ■ — ' — ' — i — ' — ' — ' — i — ' — ' — ' — i — ' — ' — ' — i — 1 T ' — i — ' — ' — ' — i — ' — ' — ' — i — ' — ' — ' — i — ' — ■ — ' — i — ' — ■ — ' — r 




7.6 7.8 8.0 8.2 8.4 8.6 7.6 7.8 8.0 8.2 8.4 8.6 

log(nf7 IvUkpc 2 ) log(n" ' M B kpc z ) 



Fig. 8. — Specific star formation rate ip* as a function of stellar mass surface density fi* for a subsample of the OPTICAL sample with 
7.8 < log^,) < 9. and M* > 10 9 ' 5 M . The left panel shows uncorrected values of , while the right shows the corrected values of ip*. 
Here again, the scatter is reduced and, notably, the slope of the relation is altered w.r.t. the uncorrected quantities. The median values of 
ip* for bins of equal size in fi* are shown as filled circles, with the bars depicting the interquartile range. The notable increase in scatter at 
high values of /i* as well as the increase in ip t may be caused by contamination from nuclear starbursts. The linear gray-scale shows the 
number density of sources at that position, with the same scale having been applied in both panels. 



in M* and /j,* according to the limits of our present 
calibration of the t b — /i* relation as defined in § 4. In 
order to establish attenuation-corrected scaling relations 
for the whole population of disk galaxies, it is crucial to 
extend knowledge of FIR/submm emission, and thereby 
calibration of the M* and /x* relation, to a representative 
subset of the entire population of rotationally supported 
disk galaxies, including dwarf galaxies. 



7. SUMMARY AND OUTLOOK 
We have presented a correlation between the face-on 
B-band dust-opacity t b and the stellar mass surface 
density ^* for normal late-type galaxies for a range in 
fi* (in units of M A:pc _2 ) of 7.6 < log(^*) < 9.0. Using 
the attenuation-inclination relation for rotationally sup- 
ported late-type galaxies we have demonstrated that the 
values of t b estimated by means of this correlation suc- 
cessfully correct statistical samples of late-type galaxies 
for dust attenuation. In order to apply the correlation 
to complete samples of spiral galaxies in terms of M*, 
the present depth of FIR/submm data used to derive 

the t b - /i» correlation means such samples need to be 
limited to M» > 10 9 5 M©. While the selection method 
employed to select spirals depends on the availability 
of NUV and optical data, only photometric optical 
data is needed to estimate the values of t b , as the 
required Sersic fits are performed on optical (r-band) 
imaging, and stellar masses have been derived using 
only optical data. As such this correlation presents a 
means to obtain attenuation corrections for samples 
of galaxies on the basis of optical photometric data 
alone, and is thus applicable to a large range of data-sets. 

We quantitatively demonstrate the efficacy of at- 
tenuation corrections using the t b - /z» relation in 
concert with the PT11 model through analysis of the 
attenuation - inclination relation in the NUV and the 
ip* vs. M* relation for a large sample of local Universe 
spiral galaxies. Both the inclination-dependent and the 
face-on components of attenuation are well predicted; 
for the latter we find that the scatter in the tA* vs. Af* 



relation is minimized for corrections within 10 percent 
of those estimated using the t b - /i» relation and the 
PT11 model, consistent with the large majority of dust 
residing in the diffuse ISM, rather than in opaque clouds. 
Overall our results are consistent with a general picture 
of spiral galaxies in which most of the submm emission 
originates from grains residing in translucent structures 
exposed to UV in the diffuse interstellar radiation field. 

We identify a roughly linear relation between the 
dust-mass and stellar mass of a galaxy for the range of 
M» probed by our sample to be the dominant driver 
for the Tg - [i % correlation. Combining this result with 
known empirical dependencies of the gas-to-stellar mass 
ratio and metallicity on stellar mass, we make a largely 
model-independent inference requiring the rapid (re- 
formation of dust grains after their destruction in the 
ISM, as has also been inferred from chemo-dynamical 
modelling of the dust cycle in the Milky- Way and exter- 
nal galaxies. Our requirement that the reformed dust is 
refractory and largely resides in the diffuse ISM (rather 
than in dense clouds) will further constrain models 

of the dust cycle in galaxies. We posit that the t b - 
/z* relation applied to statistical samples of late-type 
galaxies subdivided by, e.g, morphology/lateness and 
SFR, may provide a diagnostic tool to investigate the 
applicability to late-type galaxies in general of results 
pertaining to grain formation derived from studies of 
the local ISM in the Milky Way. 

At the moment, the determination of t b using modified 
black body functions as well as the uncertainty, intro- 
duced by the effects of diffuse dust and the bulge-to-disk 
ratio, on the measured scale lengths of galaxies are major 
contributors to the uncertainties in the results presented. 
On the other hand, however, the approach chosen in this 
paper which gives rise to these uncertainties, also ensures 
that the result given by Eq. 5 is not strongly model de- 
pendent, and is straightforward to implement in practical 
applications. In the future we will return to this topic 
using a larger data set and compare these results with 

such obtained using values of t b and s ^, r stemming 
from a fully self-consistent radiation transfer modelling 
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approach (cf. PT11). Future work will also focus on ex- 
tending the range in ^t* (hence also in M*) over which 
the relation is applicable using fully self-consistent radia- 
tion transfer modelling, two component Sersic fits (using 
future higher resolution photometry as the data becomes 
available from VST KIDS, VISTA VIKING) and a larger 
data set in the form of the entirety of H- ATLAS data. 
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APPENDIX 

APPENDIX A: THE RELATION BETWEEN DISK OPACITY, FIR/SUBMM FLUX DENSITY, AND 
DISK SCALE LENGTH IN TERMS OF THE PT11 MODEL 

In the PT11 model, the opacity of the disk of a spiral galaxy is determined by the mass distribution of the diffuse 
dust component. This is modeled as the sum of two exponential disks and has been constrained by the reproducible 
trends found in the radiation transfer analysis of the galaxy sample of Xilouris et al. (1999). For such an axisymmetric 
distribution of diffuse dust the face-on optical depth r„(r) at a given frequency v and a given radial position r is 
related to the dust surface density S(r) and the dust spectral emissivity n u as T v {r) = H{r)n v . Accordingly, for each 
disk i in the PT11 model, T v ^{r) (the face-on optical depth at frequency v and radial position r of the disk i) can be 
expressed as: 

T v ,i(r) = S ,iKrof/(^)exp ( — ) = T ,reff{v)exp I — J , (Al) 



where £o,i is the central dust surface density of the disk i, r s ,d,ref,i is the scale length of the disk i at a reference 
wavelength, K ro f is the dust emissivity at a reference frequency, and f(v) describes the frequency dependence of the 
dust emissivity given by the Weingartner & Drainc (2001) dust model {f{v) is not analytically known). Clearly, in 
this model geometry, the value of So,, is proportional to the mass of dust in the disk i and inversely proportional to 
the area of the disk, respectively the scale length squared, i.e: 



M dv 



s,d,ref,i 



(A2) 



Thus, with the opacity of the PT11 model consisting of the sum of two such exponential disks, the optical depth at a 
given wavelength and position can be fully expressed in terms of the central face-on density of dust, respectively the 
face-on opacity in a reference band (the B-band at 4430A for PT11) as: 



(A3) 



,d.B 



where K = 1.0089 pc 2 kg 1 is a constant containing the details of the geometry and the dust model of Weingartner & 
Draine (2001). 

Following the PT11 model, the total mass of dust M<i ust in a galaxy is given by 

^T^r=(i+ojc (A4) 



M dust = M** + 



where M^ st is the mass of diffusely distributed dust and M^"™ p is the mass of dust in self-shielded clumps, not 
partaking in the attenuation of optical emission. PT11 find the mass fraction of these clumps to be low (~ 10 — 15% 39 ). 
Given the uncertainties on the measurement of dust masses in comparison to the likely value of £, £ may be neglected 
to obtain: 

, (A5) 

,d,B ' s,d,B 

ducive to the formation of ices with greater emissivity coefficients, 
this estimate likely represents a upper bound 



i\/T diff 

39 This assumes the emissivity of dust in clumps is the same as 
that of diffuse dust. As cold self-shielded environments are con- 
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i.e. Eq. 1 in §2. 

Estimating from observable quantities requires several assumptions in order to re-express Eq. A3 in terms of 
observables. Under the assumption that the dust emission in the FIR, i.e. at wavelengths longwards of 100 /xm, can 
be approximated by a modified Planckian with emissivity /3, the total mass of dust in the galaxy can be expressed as: 

M M^em) S v {v oh )Dl{z)v caX 

A^dust = -j, = — — — „ — — , (AO) 



fe) B(u em ,T ) (1 + z) 1+ ^y ob m + *Kb,T ) 



where L u (u em ) is the luminosity density at the frequency v cm related to the observed frequency as v a ^ = (l + z)v t 
S v {v b) is the observed flux density, k v ^ x is the emissivity coefficient at the frequency v ca ,\, B(u, T) is the Planck function 
evaluated at frequency v and temperature T, T is the restframe temperature of the source, z is the redshift of the 
source, and Dl(z) is the source's luminosity distance. 

Similarly, the physical scale-length of the disk r S: d.B can be expressed as an angular size 9 s ,d,B as: 

r s ,d,B = e s , d , B D A (z) = e aAB f h{z) , (A7) 

I 1 + z ) 

where Da(z) is the angular diameter distance corresponding to the redshift z. Using Eqs. A4,A6, and A7 Eq. A3 can 
be expressed as: 



t 



B (1 + (l + zy+P Kvc y oh B((l + z)v oh ,To)Ol d , B Dl(z) 

K (1 + zf-V S u (v oh ) 

(1 + 0^-Kb7 2 B « l + *Kb. T o) d ld,r 

(i + zf-e s v {u oh ) 

A B((l + z)v oh ,T ) 6l ±r ' ^ j 

corresponding to Eq. 2 in §2, with the fixed geometry of the PT11 model being used to re-express r s ^ >B as f s ,d .b = 
li"a,d,r (although we set £ = in the work presented here for the purpose of determining dust masses, we have chosen 
to include the factor (1 + £) in the derivation presented here for purposes of completeness). 

Although the approximation of the dust emission from a galaxy by a single temperature modified Planckian is a 
reasonable assumption at FIR/submmm wavelengths, real galaxies will tend to have a range of components of different 
temperatures and the temperature derived will correspond to a luminosity weighted average temperature. Furthermore, 
the emissivity of the dust model of Weingartncr & Drainc (2001) is only approximately a modified Planckian with 
a fixed emissivity /3, and the actual mass fraction of dust in clumps is not known and difficult to constrain, as the 
emissivity in these regions may vary with respect to that in the diffuse medium. 

Here we have attempted to take these effects into account in first oder by empirically determining the numerical value 
of A using the radiation transfer solutions to the Xilouris et al. (1999) galaxy sample, in particular NGC891. For a 

known source with = rof , # r ,s,T-,ref = #r?s,r> ^(^ob) = Sl el (u b), z = z Te f, and T = Tq°' Eq. A8 can be used to 
identify A as: 

A _ f ^r,s,r,rcf B((l + Z rc f )v oh , T^ ci ) 

From the analysis of the Xilouris et al. (1999) galaxy sample, in particular NGC891 as presented in PT11, we 
obtain A = b.939 • 10~ 13 arcsec 2 J Jy _1 s _1 Hz -1 m~ 2 ster -1 using A ob = 250 ^m, ref = 3.5, 8 r , s ,r,rcf — 116", 
S*™V 250 ) = H5 Jy, and T rof = 20.74 K at a distance of 9.5 Mpc. 

This empirical calibration implicitly takes the mass fraction of clumps as assumed in the PT11 model into account, 
hence derived dust masses may be expected to be slightly underestimated (< 10%). 



Finally we wish to draw attention to the fact that the grain absorption cross sections in the UV and FIR of the 
Weingartner & Draine (2001) model have been empirically constrained with respect to the hydrogen gas column 
through measurements of extinction and emission of diffuse dust in the Milky Way. Thus the values of opacity are 
empirically constrained per unit hydrogen column, i.e n v — k^h- The derivation of dust masses, as e.g. given above, 
however, requires the absorption cross sections to be expressed per unit grain mass, i.e n v = K„,m- With the conversion 
unit hydrogen column to unit grain mass being relatively uncertain, the cross sections in the UV and FIR/submm 
are much more tightly constrained with respect to each other, than their absolute values. For example, as noted in 
Draine et al. (2007), the value of K v , m for the model of Weingartner & Draine (2001) requires more heavy elements 
than are expected to be available and may easily overestimate the mass of dust by a factor of ~ 1.4. In terms of the 
analysis presented here, such an overestimate will only affect absolute values such as dust masses, while leaving the 
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predicted attenuations unaffected. 



APPENDIX B: THE RELATION BETWEEN APPARENT AND INTRINSIC SIZES 

Spiral galaxies are fundamentally multi-component systems, consisting, to first order, of a disk and a bulge. 
Nevertheless, their light profiles are often fitted using single Sersic profiles, especially in the case of marginally resolved 
systems. In order to link the observed sizes, i.e the effective radius, to the intrinsic length scales of the disk and the 
bulge, multiple factors must be considered. 

While the ratio between effective radius and scale-length for a simple exponential disk is 1.68, the ratio between the 
effective radius of a single Sersic profile fit to a bulge + disk system and the scale-length of the disk component will 
decrease as the importance of the bulge increases. 

Conversely, the presence of diffuse dust in a late-type galaxy will influence the measured sizes of these objects if the 
surface density of diffuse dust possesses a radial gradient. Under these circumstances the apparent size measured 
will tend to be larger than the intrinsic size. The severity of this effect depends on both the value of t b and the 
inclination of the disk id, and is sensitive to the details of the dust geometry in the galaxy. Furthermore, as the degree 
of attenuation caused by diffuse dust varies as a function of wavelength, the effect will also be wavelength dependent. 
This effect has been quantitatively predicted for pure disk systems (Mollcnhoff, Popcscu, & Tuffs 2006; Pastrav et al. 
2013), and has been observed in the wavelength dependence of galaxy sizes (e.g., Kelvin et al. 2012; Haufiler et al. 
2012). 

A joint consideration of these effects, investigating the combined dependencies of the ratio between observed single 
Sersic effective radius and and the scale- length of the disk component on wavelength, bulgc-to-disk ratio B/D, 

inclination id, and t b has been performed by Pastrav et al. (in prep.). Fig. 9 shows the r band ratio as a function 

of inclination for 4 values of B/D at fixed t b (top) and 4 values of t b at fixed B/D (bottom). Pastrav et al., have 
performed this analysis using synthetic images of galaxies created using the same geometry assumed in PT11 and, 
accordingly, the use of these corrections is entirely consistent with the use of the PT11 radiation transfer model. 



Determining r^ s and t b 

Eq. 2 enables the determination of t b based on the observed FIR flux S250 and the angular size corresponding to 
the disk scale- length in the r band s , d ,r- Taking the corrections into account Eq. 2 can be expressed as 

T f B oc-^R(r f B ,i d ,B/D), (Bl) 

"e,ss,r 

where R{r B ,id,B / D) is the inverse of the ratio between r e , ss , r (the physical effective radius in the r band obtained 
from the single Sersic fit) and r s ,d, r as derived using the ratios of Pastrav et al. These ratios are provided in tabulated 

form, and we have interpolated the ratios in X,id, and B/D and have fit the t b dependence using a cubic spline. Using 

this spline Eq. Bl is solved numerically, obtaining the values of t b and r Sy d, r for the galaxy. 

In determining the values of t b for the OPTICAL+FIR sample we have used a value of B/D = 0.33, representative 
of the massive spirals in the sample (Graham & Worley 2008). We caution, however, that the value of B/D is a 
major source of uncertainty, which will be addressed in future work as and when higher resolution imaging, enabling 
morphological decompositions of the bulge + disk, becomes available. 
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Fig. 9. — Top: Ratio between effective radius observed fitting a single Sersic profile (r e ,ss) and the scale-length of the disk component 
(r a ,d) as a function of inclination, for four values of B/D; 0.25 (dotted crosses), 0.33 (dashed stars), 0.4 (dash-dotted diamonds), and 0.5 
(triple-dash-dottcd triangles) observed in the r band. Bottom: Ratio between effective radius observed fitting a single Sersic profile (r" a ) 
and the scale-length of the disk component (r]f) as a function of inclination, for four values of Tg at B/D = 0.33; 0.5 (dotted crosses), 1.0 
(dashed stars), 2.0 (dash-dotted diamonds), and 4.0 (triple-dash-dottcd triangles) observed in the r band. Notice the ratio of order unity 
for values of B/D corresponding to (massive) spiral galaxies, rather than 1.68 as expected for pure disk systems. Data from Pastrav et al. 
(in prep.) 
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